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' It has been suggested that an information geometric view of statistical mechanics in which a metric 

is introduced onto the space of parameters provides an interesting alternative characterisation of 
the phase structure, particularly in the case where there are two such parameters - such as the Ising 
model with inverse temperature /3 and external field h. 

In various two parameter calculable models the scalar curvature 1Z of the information metric has 
been found to diverge at the phase transition point (3 C and a plausible scaling relation postulated: 
■R. ~ |/3 -Pa] " 2 . For spin models the necessity of calculating in non-zero field has limited analytic 
consideration to ID, mean-field and Bethe lattice Ising models. In this paper we use the solution in 
field of the Ising model on an ensemble of planar random graphs (where a = —1, f3 — 1/2, 7 = 2) 
to evaluate the scaling behaviour of the scalar curvature, and find 72. ~ \f3 — fic\~ 2 ■ The apparent 
discrepancy is traced back to the effect of a negative a. 
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! I. GENERALITIES: THE INFORMATION GEOMETRY 

a 

q ■ Various authors, motivated by ideas in parametric statistics JlJ, have discussed the advantages of taking a geomet- 
rical perspective on statistical mechanics ||, [|, |J ^, J?], ||. The "distance" between two probability distributions 
in parametric statistics can be measured using a geodesic distance which is calculated from the Fisher information 
matrix for the system. In a statistical mechanical context the probability distributions of interest are Gibbs measures 
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p(x\9) = exp [-J2ff'H l (x)-\nZ{e)\ , (1) 



i=l 



CN 
> 
CO 

in 

■ where the x characterise the state of the system (e.g. spins), the Hi{x) are the various terms in the Hamiltonian, the 
Z(6) is the normalising partition function and the l are the various parameters such as the inverse temperature (3, 
. the external field h, etc. 

The manifold A4 of parameters is endowed with a natural Riemannian metric, the Fisher- Rao metric Q, which 
measures the distance between different configurations. For a spin model in field A4 is a two-dimensional manifold 
parameterised by (6' 1 ,6 |2 ) = (f3,h). The components of the Fisher-Rao metric take the simple form Gtj — didjf in 
this case, where / is the reduced free energy per site and di = d/d9 l . A natural object to consider in any geometrical 
approach is the scalar or Gaussian curvature which may be calculated as 
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d$f dp h f dpdif , ( 2) 

dphf dpdlf dlf 
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where G = det(Gy). 

It is worth remarking that, unlike most standard statistical mechanical observables, the curvature TZ depends on 
third order derivatives. Nonetheless, a plausible scaling relation has been advanced for TZ in the critical region. 
The hypothesis, on dimensional grounds, is that the curvature depends on the correlation volume for a second-order 
transition TZ ~ £ d , where £ is the correlation length and d is the dimension of the system. This is reasonable since £ 
is the only physical scale in the system near criticality. 

Combined with hyperscaling, vd = 2 — a, and standard scaling assumptions this leads to 



(3) 
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In the above a is the standard exponent characterising the scaling of the specific heat, so consideration of TZ clearly 
offers a way of determining critical exponents in a non-standard manner. 

Analytic determination of 1Z in spin models has been limited by the necessity of carrying out calculations in field. 
One case where this is possible is the ID Ising model 0, where the curvature was calculated to be 

1Z = 1 + rr 1 cosh/i (4) 

with r\ = \J sinh 2 h + e~ 4 ^. The ID Ising model can be thought of as having a zero-temperature transition, so looking 
at h — 0, (3 — * oo we see that 1Z ~ e 2/5 , corresponding to the expected a = 1. Similarly, it is possible to obtain 
an expression for the scalar curvature for the Ising model on a Bethe lattice ||, where the scaling behaviour is also 
verified. Both of these examples have unsatisfactory aspects - the ID Ising model has no real transition and the Bethe 
lattice Ising model is mean field in nature. 

Given the relative paucity of models which are soluble in field any further explicit calculations would be welcome, 
particularly in a non-mean-field model with a genuine finite-temperature phase transition. In the sequel we discuss 
one such case, the Ising model on dynamical planar random graphs. 



II. PARTICULARS: THE ISING MODEL ON PLANAR GRAPHS 



The solution of the Ising model on an ensemble of $ 4 (4-regular) or <E> 3 (3-regular) planar random graphs was 
first presented by Boulatov and Kazakov ||, 10 1, who were motivated by string-theoretic considerations, since the 
continuum limit of the theory represents matter coupled to 2D quantum gravity. They considered the partition 
function for the Ising model on a single n vertex planar graph with connectivity matrix G"j 



then summed it over all n vertex graphs {G™} resulting in 

\ c (G n ,P,h) , 

before finally forming the grand-canonical sum over graphs with different numbers n of vertices 
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where c = exp(— 2/3). This last expression could be calculated exactly as matrix integral over N x N Hermitian 
matrices, 



Z = - log / T>4>\ V(j> 2 exp -Tr 



l(<fi + -cfafa +e- h <&) 



(8) 



where the N — > oo limit is to be taken to pick out the planar diagrams and the potential appropriate for <E> 4 (4-regular) 
random graphs has been shown. 

When the matrix integral is carried out the solution is given in parametric form by 
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where the function g(z) is 



q(z) = -c 2 z 3 + - 
yy ' 9 3 



dt 
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and B = 2[cosh(ft) - 1]. 

In the thermodynamic limit the free energy per site is given by 
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where zq = Zq{[}, h) is the appropriate low- or high-temperature solution of g'(z) — 0. When h = this may be solved 
in closed form, and although the solution is not available explicitly for non-zero h it can still be developed as a power 
series in h around the zero-field solutions in order to obtain expressions for quantities such as the energy, specific 
heat, magnetization and susceptibility. It was found that the critical exponents were given by a = —1, (3 = 1/2, 
7 = 2, so the transition was third order with, intriguingly, the same exponents as the 3D spherical model on a regular 
lattice 0. 

If we carry out a perturbative expansion for the high-temperature solution, which is symmetric in h and hence a 
series in even powers, we find 

1 (w- l)(2u 2 - 2u+ 1) 2 

z - 1 77 Tyi "> 



u (2u-l 



(u - l)(2u 2 - 2u + l)(4u 5 - 10?/ + 10u 3 - 5m 2 + 5w+1) 4 
+ 24(2u-l)9 + '" ' ( } 

where the coefficients in the series are most naturally expressed in terms of u = exp(— 0) = t/c, as above. 

III. GENERALITIES: SCALING OF THE SCALAR CURVATURE 

The expected scaling form of the various components of 1Z for a generic spin model in field is discussed at some 
length in Q] , and we now recapitulate these results briefly for comparison with the specific results for the Ising model 
on planar random graphs in the next section. The starting point is the scaling form of the free energy per spin near 
the critical point, 

f(e,h) = \- 1 f(e\ a ',hX^) , (13) 

where e = f3 c — (3 and a e ,a,h are the scaling dimensions for the energy and spin operators. For e > 0, i.e., in the 
unbroken high-temperature phase, we can use standard scaling assumptions to write this as 

f(e,h)=e 1 / a ^+(he- a ^), (14) 

where ip + is a scaling function and we also define A — l/a e and C — —ah/a e for later convenience. In terms of the 
standard critical exponents A = 2 — a and A + C = (3. 

This generic scaling form can now be substituted into equ. (^|) to find the scaling behaviour of the various components 
and the scalar curvature (0) itself near criticality (i.e. h = 0, e — ► 0), 



A(A- l)e A - 2 V;+(0) e A+2C V+(0) 

—A(A — I) (A — 2)e A ~ 3 ip + (0) -{A + 2C)e A+2C - 1 ^l(0) 

~{A + 2C)e A+2C - 1 ^l(0) e A+3C ^{0) 



(15) 



where the scaling of the metric determinant is 

G = A(A - l)e 2A+2C ^ + (0)i/C(0) • ( 16 ) 
Expanding the determinant one finds two terms of similar order contributing to give 

„ (A + 2C)[(A + 2C)-(A-2)} A 

2A(A-l)V> + (0) 6 [U) 

or, translating back to the standard critical exponents, 

2(2-a)(l-a>M0) e ' 

The discussion in pj was intended to be as general as possible, one should note that for Ising-like models with a ±h 
symmetry all odd derivatives of the scaling function w.r.t. h will vanish so 9|/ = rather than e A+3C ^p + (0). This 
does not affect the stated scaling relations. 

However, one feature of these scaling relations does have an impact on the observed scaling for the Ising model. 
Generically one expects that 91/ = A(A — l)e A ~ 2 ip + (0), which contributes to both the metric and the determinant 
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involved in calculating TZ. If A > 2, i.e. a < 0, this naively suggests a vanishing 92/ at criticality, which will in 
general not be the case. There would instead be a contribution from a regular term, which would give a constant at 
the critical point. Having such a constant term modifies the scaling form of TZ in the case a < 0, A > 2 to 
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2G 2 



A(A- 1)0(0) 
-A(A-l)(A-2)e A -V + (0) 






C A+2C-1, 



C A+2C 



-{A + 2C)e 



vC(o) 

A+2C-1 



<Mo) 



-(A + 2C)e A + 2G '- 1 V'+(0) e A+3U ^ (0) 

where we have denoted the constant by A(A — 1)0(0). The scaling for G is also modified to 

G = A(A-l)e A+2C (/)(0)ipl(0) . 



(19) 
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When expanded the expression for TZ contains two term which now have differing orders in e. The leading term for 
A > 2, the case which we are interested in, is 

(A + 2C) 2 _ 2 
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2A{A- 1)0(0) 



(21) 
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2(2-a)(l-a)0(O) 



(22) 



so the critical exponent a no longer appears in the scaling exponent. 

By virtue of the Boulatov and Kazakov solution the Ising model on planar random graphs allows us to explicitly 
confirm these observations, as we see in the next section. Since a = — 1,/3 = 1/2,7 = 2, we have A = 3, C = — 5/2 
and the modified discussion of scaling should apply. 



IV. PARTICULARS: THE SCALAR CURVATURE FOR ISING 



We can now take the series expansion for zq from equ. (|12]), insert this into g(z) and use the resulting expression 
for / in equ. ( pi] ) to calculate the various terms that appear in the scalar curvature TZ as power series in h 2 . We find 
that the leading terms at h = 0, with e u = u — u cr — e/2 + . . . and u cr = 1/2, and using /?, h as co-ordinates are 



TZ = - 
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2G2 
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675 
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(23) 



The determinant of the metric scales as G - 

TZ 



88 -2 
375 fc " 

225 n 

704 u 



so the final scaling expression for the scalar curvature is 
225 , 



176 



(24) 



A glance back at equs. (|T^)-(|2^) shows that the modified scalin g fo r A > 2 that these incorporate is, indeed, followed 
for the individual components in equ. (|l9|), the metric in equ. (pfj|) and the scalar curvature itself in equs. (21 ) and 



It is an easy matter to calculate TZ for any u when h is small, using the expansion for zq in equ. (fL2J). Writing 



TZ = TZ + TZ 2 h 2 



(25) 



the first two coefficients are given by 

(6m 8 



TZo = 



43m 7 + 357m 6 + 1265m 5 + 2123u 4 + 1841u 3 



783u 2 + 75m + 3) 



2(6m 5 + 27u 4 + 56u 3 + 54m 2 + 18m + 3) 2 (2m - l) 2 
x (2m 2 + 2m + 1)(m+1)(m 2 + 1) 



(26) 



and 
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TZ 2 = -(u+l)(u 2 + l)u{u + 2) [144m 18 + 1008? 



17 - 3276m 16 - 31180m 15 



79106m 1 



- 129786m 13 - 135424m 12 - 92093m 11 - 78645m 10 
+ 245658m 7 + 410788m 6 + 328760m 5 + 139986m 4 + 33183m 3 
+ 45]/ [(6m 5 + 27?i 4 + 56m 3 + 54m 2 + 18m + 3) 3 (2m- l) 7 ] 



37499m j + 54941m b 



5331m 2 + 765m 



(27) 
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which give the scaling of equ. (|24|) when u is set equal to 1/2 + e u . 

In Fig. 1 we have plotted TZ close to u cr — 1/2 using a series correct up to 0(h e ) terms. The scaling region in 
h is very narrow, with the approximation to TZ rapidly giving large negative values outside this region due to the 
increasingly strong divergences in the series coefficients as u — > u cr for increasing order in h. This turnover is just 
visible on the edges of the plotted surface. The sensitivity to h would have to be carefully handled in any numerical 
investigations of TZ. Within the domain of validity of the expansion in h it appears that TZ is positive. It has been 
remarked || that 1Z is positive in the thermodynamic limit for Ising models when the parameters take physical 
values, with the only divergence being at the critical point, and the Ising model here provides another example. This 
feature is apparently not universal, calculations of TZ for the one-dimensional Potts model |l2| and field theories || 
do not give positive curvatures throughout the physical parameter space. 

It has been observed that the line h = is a geodesic of the metric for the one-dimensional Ising and Potts models 
p2| . The geodesic equations using co-ordinates 0, h are given in general by 

^ + r% f ,V f, V f> + 2T^ h V?V h + Ti h V h V h = \{a)V p , (28) 
dV h 

— + r^W" + 2Y h 0h V p V h + T h hh V h V h = \{s)V h , (29) 

where s parameterises the flow lines, V 13 = d/3/ds, V h = dh/ds, the T are the Christoffel symbols and A(s) allows for 
the possibility of a non-affine parameter choice. 

A vector field with a flow line along h = has V h — 0, so in this case equs. (^8|) and ( p9| ) reduce to 

■ + { T "p)\ h =o VPvP = ( 3 °) 
^)L^ = . (31) 



ds 



The first of these equations for /3(s) always has a solution and the second requires ( Tgg )L_ n = 0. This is satisfied for 



0PJ\h=O 

the Ising model on planar random graphs because the Christoffel symbol vanishes at h = for the same reason as 
in the one-dimensional models - the off-diagonal components of the metric dpd^f are O(h) in all cases. We therefore 
find that for the Ising model on planar random graphs h = is also a geodesic line in the /?, h plane. 

We close with a remark on the infinite temperature (T— >oo, /3 — > 0, u — > 1) limit of TZ when h = 0. In this limit 
TZ was found to be 2 for the one-dimensional Ising model || and z/2 for the Ising model on a z co-ordinated Bcthc 
lattice ||. Here we find that TZ(T = oo) = Yg|y = 2.415...., so if we accept the suggestion in |jj that TZ — TZ(T = oo) 
should be taken as a measure of fluctuations caused by the spin interactions the correct measure of the deviation from 
ideal paramagnetism for the Ising model on planar random graphs is TZ — -f§f"[- 



V. CONCLUSIONS 



We have calculated the scaling behaviour of the scalar curvature of the Fisher-Rao metric for the Ising model on 
planar random graphs using the exact solution of |l0| combined with a perturbative expansion in the external field 
h. Although a = —1 for this model, TZ ~ e~ 2 rather than the TZ ~ e~ 3 one might have expected naively from general 
scaling arguments. This discrepancy was traced back to the effect that a negative value of a had on the scaling of the 
various components of the metric and the terms which contributed to TZ. 

Various qualitative features of the calculated TZ tally with earlier observations of one-dimensional and mean-field 
Ising models. It is positive (within the domain of applicability of our semi-perturbative calculation) and diverges 
only at the critical point. The zero-field line is seen to be a geodesic, just as for the one-dimensional Ising and Potts 
models. 

It would be an interesting exercise to calculate TZ for other models where some form of solution in field was accessible. 
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FIG. 1: A plot of TZ close to u cr = 1/2. Note that the external field h is scaled by a factor of 10 5 and so covers a very narrow 
range. 



